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.£h , 1 INTRODUCTION 



ABSTRACT 

We use the statistics of strong gravitational lensing from the CLASS survey to impose 
constraints on the velocity dispersion and density profile of elliptical galaxies. This 
approach differs from much recent work, where the luminosity function, velocity dis- 
persion and density profile were typically assumed in order to constrain cosmological 
parameters. It is indeed remarkable that observational cosmology has reached the point 
where we can consider using cosmology to constrain astrophysics, rather than vice 
versa. We use two different observables to obtain our constraints (total optical depth 
and angular distributions of lensing events) . In spite of the relatively poor statistics and 
the uncertain identification of lenses in the survey, we obtain interesting constraints 
on the velocity dispersion and density profiles of elliptical galaxies. For example, as- 
suming the SIS density profile and marginalizing over other relevant parameters, we 
find 168km/s < cr* < 200km/s (68% CL), and 158km/s < cr* < 220km/s (95% CL). 
Furthermore, if we instead assume a generalized NFW density profile and marginalize 
over other parameters, the slope of the profile is constrained to be 1.50 < j3 < 2.00 
(95% CL). We also constrain the concentration parameter as a function of the density 
profile slope in these models. These results are essentially independent of the exact 
knowledge of cosmology. We briefly discuss the possible impact on these constraints 
of allowing the galaxy luminosity function to evolve with redshift, and also possible 
useful future directions for exploration. 



3 



The statistics of strong gravitational lensing has repeatedly 
been advertised and used as a probe of cosmology (e.g., 
Turner, Ostriker & Gott 1984, Hinshaw and Krauss 1987, 
Fukugita et al. 1992, Krauss and White 1992, Kochanek 
1995, 1996, Cooray, Quashnock & Miller 1999, Chiba & 
Yoshii 1999, Cheng & Krauss 1999). The sensitivity of lens- 
ing counts to Q,m and Qa, the energy densities in matter 
and the vacuum component relative to the critical, comes 
mostly from a volume effect: higher Ha implies bigger co- 
moving volume for a fixed redshift, leading to the higher op- 
tical depth for lensing. Using knowledge about the luminos- 
ity function of galaxies and their density profiles, many au- 
thors have used lensing statistics to constrain cosmological 
parameters. For example, Fukugita and Turner (1991) first 
constrained the vacuum energy density to be less than about 
90% of the critical energy density (Qa < 0.9) at 95% con- 
fidence level (hereafter CL). Subsequently this was followed 
by Kochanek (1995, 1996), who claimed an upper limit on 
the vacuum energy density (Qa < 0.66 at 95% CL). Krauss 
and White (1992) and later Chiba and Yoshii (1999) and 
Cheng and Krauss (1999) used a different choice of galaxy 



parameters and demonstrated that a flat vacuum-energy- 
dominated universe could be favored. Similar analyses have 
been performed by Im et al. (1997), Cooray, Quashnock & 
Miller (1999), Waga & Miceli (1999), and all typically favor 
the ACDM cosmology. Cheng and Krauss (1999, 2001) also 
explored how uncertainties in the choice of galaxy parame- 
ters could result in vastly different constraints on cosmology, 
although they argued for a choice that ultimately favored 
a flat, vacuum-energy-dominated cosmology. It has also re- 
cently been argued that strong-lensing statistics from ongo- 
ing surveys like the Sloan Digital Sky Survey (SDSS) might 
impose interesting constraints on the equation-of-state ratio 
of dark energy w (Cooray & Huterer 1999); constraints on w 
from lensing have already been claimed by Sarbu, Rusin & 
Ma (2001) who used the statistics of the JVAS/CLASS sur- 
vey to obtain w < —0.4. Similar results have been obtained 
very recently by Chae et al. (2002). 

Given the notoriously poor statistics of strong lensing 
surveys thus far — the total number of gravitational lenses 
is of order fifty, and the largest homogeneous survey (which 
we use in this work), JVAS/CLASS, currently has a total 
of only 17 events — combined with the existing galactic lu- 
minosity function uncertainties, it is not clear how seriously 
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one should take any constraints on cosmology derived from 
strong lensing statistics. To constrain cosmological parame- 
ters using lensing statistics one has to deal with the strong 
dependence of the results on the lens profile, the density 
dispersion of galactic dark matter, the number density of 
galaxies as a function of redshift, and observational effects 
due to magnification bias and the selection function of the 
survey. 

In this work, we exploit this sensitivity to reverse the 
traditional methodology. Since lensing statistics are, on the 
whole, much more sensitive to astrophysical than cosmolog- 
ical parameters, we wish to utilize existing surveys to probe 
the properties of lensing galaxies rather than cosmology. We 
are aided in this effort at this time because independent 
probes of cosmological parameters have recently converged 
rather tightly on a single cosmological model: a flat dark 
energy dominated universe with Ode ~ 0.7, and Q.m ~ 0.3. 
As these parameters currently seem to be more tightly con- 
strained that the galaxy parameters described above, now 
seems an opportune time to use cosmology to constrain as- 
trophysics, rather than vice versa! 

Some efforts along these lines have already been ex- 
plored, as new and better lensing data, especially the 
JVAS/CLASS survey, have appeared. In particular, several 
investigations have been undertaken to constrain the nature 
of galaxy clustering in the CDM paradigm. Keeton (2001) 
used the statistics of JVAS/CLASS lenses to indicate that 
CDM galaxies are too concentrated to agree with the lensing 
statistics, while Keeton & Madau (2001) used the absence 
of wide-separation lenses in the CLASS survey to impose 
an upper bound on the concentration of dark matter halos. 
Takahashi & Chiba (2001) consider lensing by both singular 
isothermal sphere (SIS) and Navarro-Frenk- White (NFW) 
profile galaxies, and find that the lack of observed large- 
angle separation lenses indicates that the density profile is 
not too steep (/3 < 1.5, with p(r) oc r~^). Oguri, Taruya 
& Suto (2001) obtained a similar result by using the statis- 
tics of tangential and radial arcs. Conversely, Rusin & Ma 

(2001) use the absence of detectable odd images to set a con- 
straint on the surface density of lensing galaxies, and con- 
clude that lenses cannot have profiles much shallower than 
an SIS {/3 > 1.8). Wyithe et al. (2001) and Li & Ostriker 

(2002) considered lensing by objects with both SIS and gen- 
eralized NFW (CNF W) density profiles. They computed op- 
tical depths, image separations and magnification biases. In 
particular, Li & Ostriker, extending the earlier work of Kee- 
ton (1998) and Porciani & Madau (2000) argued that, in or- 
der to explain the large number of observed small-separation 
lenses and the lack of large-separation events (compared to 
predicted distributions for lensing by clusters), the favored 
galaxy cluster profile seems to be the combination of SIS 
(when M < 1O 13 M ) and NFW (when M > 1O 13 M ). 

Here we carry out a related analysis, with the aim of 
constraining the nature of individual galaxies rather than 
clusters. For this purpose we shall assume the "concordance" 
values for the cosmological parameters (e.g. Krauss 2000): 
Am = 1 — !Jde = 0.3, w = — 1 and h = 0.7, where Qm and 
JIde are energy densities in matter and dark energy relative 
to critical, w is the equation of state ratio of dark energy, 
and Hq — 100 h km/sec/Mpc. We will show that our results 
are extremely weakly dependent on the assumed cosmology 
(in particular, knowledge of Qm). 



2 THE DATA 

Although more than 60 multiply imaged quasars and radio 
sources are known, they come from different observations 
with different sensitivities and selection functions, which 
makes an accurate computation of the expected number of 
lenses very difficult. Therefore, it is imperative to have data 
from a single well-understood survey with information on 
the source population. In this work we use the most com- 
plete homogeneous sample of lenses provided by the Cosmic 
Lens All-sky Survey (CLASS; Myers et al. 2002, Browne et 
al. 2002) , which extended the earlier Jodrell-Very Large Ar- 
ray Astrometric Survey (JVAS; Patnaik et al. 1992a, King 
et al. 1999). CLASS is using the Very Large Array to im- 
age radio sources with the flux density of between 30 and 
200m Jy; candidate lensing events are followed up by Multi- 
Element Radio-Linked Interferometer Network (MERLIN) 
and the NRAO Very Large Baseline Array (VLBI). So far a 
total of about 16,000 sources were imaged by JVAS/CLASS, 
with 22 confirmed lensing events. Of these, a subset of 8958 
sources with 13 lenses forms a well-defined subsample suit- 
able for statistical analysis (Browne et al. 2002), and we use 
this subsample in our work. Table 1, essentially identical to 
Table 3 in Browne et al. (2002), shows the lenses from the 
statistically controlled subsample. We have added informa- 
tion about the identity of the lens, in particular whether it 
is a spiral galaxy, an elliptical, or formed by more than one 
galaxy (Chae 2002). 

It is well known that elliptical galaxies dominate the 
optical depth for strong lensing by individual galaxies (e.g. 
Kochanek 1993b), and as a result we concentrate on con- 
straining their parameters here. This effort is somewhat 
complicated by the fact that only six of the CLASS lenses 
are clearly identified as ellipticals and one as a spiral, while 
in other cases the identity of the lens is uncertain; see Table 
1. Furthermore, three events are due to more than one lens 
galaxy. It is crucial to choose a subset of CLASS lenses that 
includes elliptical galaxies only. It is clear that the num- 
ber of ellipticals is between 6 and 12, and that confirmed 
ellipticals outnumber spirals in ratio 6:1. The most likely 
value of the number of ellipticals is therefore somewhere 
near 11. To compute the measured number of lenses, we 
chose marginalize over the range between 6 and 12, with the 
gaussian weighting centered at 11 and variance of 5. How- 
ever, as we later discuss, the results are extremely insensitive 
to the exact choice of weighting; the reason is that the statis- 
tics are much more sensitive to the parameters we wish to 
constrain, the velocity dispersion and density profile of ellip- 
tical galaxies. For the angular separation test, we use only 
the four single* elliptical lenses (B0712+472, B1422+231, 
B1933+503, B2319+051). To test the robustness of this test, 
we alternatively assume that all unidentified galaxies are el- 
lipticals as well, and use a total of 9 single non-spiral lenses 
(the four above, plus B0445+123, B0631+519, B0850+054, 
B1152+199 and B2045+265). As discussed later on, our re- 
sults are insensitive to the exact choice of this subset. 

We wish to utilize three different observables to obtain 

* Multiple lens deflectors obviously make different predictions 
from single galaxies, and produce larger angular separations of 
images. We are interested in splittings due to single elliptical 
galaxies only. 
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Survey 


Lens 
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Zs 


/) 

U 


TT~1 
1JJ 


Reference 


JVAS 


B0218+357 


0.68 


0.96 


0.33 


S 


Patnaik et al. 1993 


CLASS 


B0445+123 


0.56 





1.33 


? 


Argo et al. 2002 


CLASS 


B0631+519 








1.16 


? 


Browne et al. 2002 


CLASS 


B0712+472 


0.41 


1.34 


1.27 


e 


Jackson et al. 1998 


CLASS 


B0850+054 


0.59 





0.68 


? 


Biggs et al. 2002 


CLASS 


B1152+199 


0.44 


1.01 


1.56 


? 


Myers et al. 1999 


CLASS 


B1359+154 


— 


3.21 


1.65 


?, m 


Myers et al. 1999 


JVAS 


B1422+231 


0.34 


3.62 


1.28 


e 


Patnaik et al. 1992b 


CLASS 


B1608+656 


0.64 


1.39 


2.08 


e, m 


Myers et al. 1995 


CLASS 


B1933+503 


0.76 


2.62 


1.17 


e 


Sykcs et al. 1998 


CLASS 


B2045+265 


0.87 


1.28 


1.86 


? 


Fassnacht et al. 1999 


JVAS 


B2114+022 


0.32/0.59 




2.57 


e, m 


Augusto et al. 2001 


CLASS 


B2319+051 


0.62/0.59 




1.36 


e 


Rusin et al. 2001 



Table 1. Thirteen lensing events from the "CLASS statistical sample" of 8958 objects (adopted from Browne et al. 2002; see also Chac 
2002). "ID" stands for identification of the lens - whether it is a spiral galaxy (s), an elliptical (c) or unknown (?); three lenses consist 
of multiple galaxies (m). 



our constraints: the overall optical depth r to a source at 
redshift z s , the differential optical depth as a function of 
angular separation, and the differential optical depth as a 
function of lens redshift. Unfortunately, the last of these 
tests is uncertain due to possible incompleteness of the sur- 
vey: higher-redshift lenses are more difficult to measure due 
to their lower fluxes; while the source redshifts are more eas- 
ily measurable for objects very far away (mainly quasars) 
and very close (mainly galaxies), and not ones at interme- 
diate distances. Because of these uncertainties, and because 
the redshift test does not add much to our constraints, wc 
decide not to use the redshift-distribution test^ 

We are therefore left with two tests, the total optical 
depth (r-test) and angular separation (dr / d9 -test) . The for- 
mer test gives stronger constraints in both SIS and GNFW 
cases. The latter test, in the SIS case, is independent of z s as 
long as z a > 0.2; henceforth, the knowledge of z s is not nec- 
essary and all single ellipticals (chosen as explained above) 
can be used for this test. In the GNFW case the knowledge 
of z 3 is required for this test, and when it is not available we 
use the mean redshift of the measured sources, z s = 2. (We 
have checked that the results change negligibly if, instead 
of z s = 2, we use the histogram of the source distribution 
from Marlow et al. (2000), which is centered at z s = 1.27 and 
has long tails.) Finally, we use the maximum lens separation 
#max as an estimator of the angular separation 9. Although 
this estimator has been widely used in the literature due to 
the fact that # ma x are readily available, we warn that the 
angle corresponding to the average image radius fitted to a 
lens model, for example, would be a better estimator. Nev- 
ertheless, we do not expect that using # max will significantly 
bias the results, given the limited current statistics. More- 
over, since higher <r* roughly corresponds to larger angular 
separations, our results may only be biased to higher a,, 



strengthening our conclusion that this parameter is smaller 
than previously quoted in the literature. 

In order to compute the expected optical depth for 
any given model, it is crucial to know the redshifts of 
source quasars and galaxies. The redshift distribution of 
JVAS/CLASS source objects has been discussed in Marlow 
et al. (2000), who spectroscopically followed up 42 sources 
at William Herschel Telescope. Most of these sources are 
quasars; with a significant admixture of galaxies at z < 1. 
The mean redshift of this subsample is (z s ) = 1.27 with an 
rms spread of 0.95. In this work we use the full histogram 
distribution of the observed subsample of sources (Fig. 2 
in Marlow et al. 2000), and assume that the redshift dis- 
tribution of the subsample gives a good representation of 
the overall redshift distribution. One has to be cautious, 
however, since the lensed sources come from a fainter pop- 
ulation than the ones in Marlow et al. (2000), and may be 
at different redshifts. The validity of this assumption has 
been examined by Chae (2002) , who reviews existing obser- 
vations and finds that the redshift distribution is expected 
not to change much at lower flux densities, corresponding to 
lensed sources. 

Finally, we will need to know a few other details re- 
garding the CLASS sample. The survey is complete at image 
separations 0.3" < 9 < 15" (Helbig 2000, Myers et al. 2002). 
All confirmed JVAS/CLASS lenses have image separations 
9 < 3". The distribution of sources as a function of the total 
flux density S is well described by the power law 



where dn is the number of sources observed in flux density 
interval dS. For JVAS/CLASS, 7} ~ 2.1 (Rusin & Tegmark 
2001). 



t Nevertheless, we have checked that the results of the redshift- 
distribution test agree with those of the other two tests. Further- 
more, for the SIS case the quantity (1 / 't)(cLt / 'dzi) is independent 
of galaxy parameters, and we used it to check that constraint 
on fijj and w is consistent with the adopted cosmological model 
f!jw = 0.3, w = —1. 



3 DENSITY PROFILE 

There is good evidence that the density profiles of dark 
halos on cluster scales depend on the halo mass (Keeton 
1998, Wyithe et al. 2001, Li and Ostriker 2002). For the 
less massive halos (M < 10 13 Mq), SIS profiles are found 
to be adequate, while for large-mass halos (M > 10 13 Mq) 



© 2002 RAS, MNRAS 000, 000-000 



4 Davis, Huterer & Krauss 



NFW profiles provides a good fit. This result is also ex- 
pected from semianalytic models, which show that objects 
smaller than M c w 10 13 Mq are subject to baryonic cooling, 
whereby baryons collapse to the center thereby enormously 
increasing the central density and lensing cross-section, and 
converting the shallow NFW profiles into the steep SIS (Rix 
et al. 1997, Kochanek & White 2001, Keeton 2001). 

Galaxy clusters tend to lead to large lens separations 
and/or extended arcs and arclets. Since 



1.271 



D s 



M 



© 



2/3 



/ /Ocrit 
V Pcrit,0 



1/3 



(2) 



(Li & Ostriker 2002) where all quantities except p C rit,o are 
evaluated at z = zi, we see that M > 10 Mq corresponds 
to 9 > 3". If we are interested in primarily lensing by indi- 
vidual galaxies rather than clusters, we should concentrate 
on image separations substantially smaller than this value. 
In CLASS, all lensing events have separations smaller than 
3". We therefore conclude that the CLASS lenses are due to 
individual galaxies and not clusters. 

Furthermore, we assume smooth, spherically symmetric 
density profiles. This assumption is widely used, and sup- 
ported by the findings that the subclumps do not greatly 
affect the total optical depth for lensing (Flores, Mailer & 
Primack 1996) and that asphericity of density profiles affects 
mostly the ratio of quads to doubles and not the optical 
depth (Rusin & Tegmark 2001). 

Our goal in this paper is twofold. First, we would like 
to constrain the galaxy velocity dispersion assuming the 
SIS profile. The SIS profile has repeatedly been used in 
the past to constrain cosmological parameters, assuming the 
Schechter function parameters and the galaxy velocity dis- 
persion to be known. We would like to reverse this process 
and see whether the previously-used <r* is still favored now 
that we have good knowledge of cosmological parameters. 

Second, we would like to constrain the density profile 
of elliptical galaxies. As argued above, only the inner parts 
of lens galaxies (a few tens of kiloparsecs from the center) 
are responsible for CLASS events. Moreover, as discussed 
in Sec. 7, there is good evidence that cores of galaxies are 
small and can safely be ignored. Therefore, it seems justified 
to adopt p(r) oc r _/3 and try to constrain f3. We do this via 
the generalized Navarro-Frenk- White profile, as described in 
Sec. 7. 



4 MODELING THE LENS: SINGULAR 
ISOTHERMAL SPHERE (SIS) PROFILE 

4.1 Number density of lenses 

Since we are interested in elliptical galaxies, we adopt the 
Schechter luminosity function (Schechter 1976) which has 

repeatedly been shown to be a good fit to the measurements^ 



t In order to compute optical depths for generalized dark matter 
distributions on cluster scales, many authors have assumed the 
Press-Schechter mass function (Press and Schechter 1976). Since 
we are interested in constraining observational properties of ellip- 
tical galaxies, and since the lens identification from the CLASS 
survey indicates that most lenses are due to individual galax- 



dL 



(L) dL = 4>. t 



exp(— L/ L t ) 



dL 



(3) 



There has been much discussion as to what values of (f>„ and 
a best describe the actual luminosity function. Typically, it 
is argued that ^>„,tot = 1-4 x 10 h Mpc for all galaxies, 
of which « 30% are ellipticals (Postman & Geller, 1984), so 
that <f)t mp = 0.6 x 10~ 2 ft 3 Mpc~ 3 ; further, a « -1 with fairly 
large uncertainties. Recently the SDSS (Blanton et al. 2001) 
claimed a more accurate determination of the local (z < 0.2) 
luminosity function; a — — 1.20±0.03 and <^>*,tot = (1.46± 
0.12) x 10" 2 /i 3 Mpc" 3 . 

To relate the luminosities to velocity dispersions, we use 
the Faber- Jackson relation (Faber & Jackson 1976) 



(4) 



where it is typically assumed that 7 w 4 for the SIS profile. 
Our principal goal is to determine the parameters 0*, a, 7 
and a*. 



4.2 Optical depth 

The optical depth for a lens at redshift zi due to a particular 
source at z s is given by 



f" a dDl n , ^3 

/ dzi-^— (1 + 2;) X 

Jo 



dzi 



L 



dL-£(L) <xsis(L, z u z.)B{L, z h z s ) (5) 

where <f> is the comoving number density of lenses, L is their 
luminosity, asis(L, zi, z s ) is their cross-section for lensing, 
and B(L, zi, z s ) is the magnification bias, describing the fact 
that lensed galaxies will be magnified, and therefore seen 
more easily, and therefore are enhanced in any flux limited 
survey. In Eq. (5) we have allowed for a general redshift and 
luminosity dependence of the number density, cross-section 
and magnification. For redshift-independent (as we first as- 
sume) 0*, a and 7, d<j)/dL depends only on L. Similarly, 
assuming that CLASS lenses are described by an SIS pro- 
file and the radio luminosity function is a power law, B is 
simply a constant (see below). 

The density profile for the SIS is given by 



p(r) 



2nGr 2 



(6) 



where a is the velocity dispersion of the galaxy. This dis- 
tribution produces an image separation of 29 e, where the 
Einstein radius 9b = 47t(ct /c) 2 Di s / D s , and Di s and D s are 
the angular diameter distances between the lens and source, 
and observer and source respectively. The cross-section for 
lensing is therefore 



o"srs = tt(9 e Di) 2 = 16n 3 



(!)' 



DiDu 
D s 



(7) 



We only consider angular separations greater than some 
minimum value # m in, since the resolution limit of CLASS is 
#min = 0.3", and multiple images with smaller separation 



ies, for our purposes the Schechter luminosity function is more 
relevant. 
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Figure 1. Dependence of the observables on the velocity dispersion cr*, assuming that all other parameters take their fiducial values. 
Left panel: The dependence of r (the shaded region is the measured value from CLASS, assuming the number of ellipticals to be 
between 6 and 12). Right panel: The dependence of (l/r)(dr/d6) (vertical lines denote measurements from CLASS, solid lines denote 
confirmed ellipticals, while dashed lines denote galaxies whose type has not been identified). The other Schcchtcr-function and cosmological 
parameters were fixed to their fiducial values from Sec. 4.3. Note that the CLASS survey is complete for 8 > 0.3"; therefore, all predicted 
quantities, such as dr/dd in the right panel, were compared to measurements only for 6 > 0.3". The quantity dr/dzi is very weakly 
dependent on <r» (and other Schechter function parameters) and is not shown. 



angles than # m i n will not be resolved. The correspondence 
between the luminosity and angular separation for an SIS 
lens is 



L = 



W a c 2 



8tt D iB a; 



7 /2 



/. 



(8) 



where c is the speed of light, so that m i n corresponds to 
SOme 1/min as the lower limit of integration in Eq. (5). 

We also need to compute the magnification bias. It is 
given by 



B = 



dn S 
dS m 



P(/i)/i 1 dfi 



dn 
dS 



(9) 



where dn/dS is the source luminosity function and P(fi) the 
distribution of total magnifications. For the power-law lumi- 
nosity function of CLASS (cf. Eq. (1)) and the distributions 
of magnifications for SIS lenses (P(p) = 8/i -3 ), the bias 
simplifies to (Sarbu, Rusin & Ma 2001) 



B(L,zi,z s ) = 4.76. 



(10) 



Finally, we will be interested in the quantity dr/dO. This 
is given by 



dr 
dO 



dr 
dh 



x — 

L(e t ) de 

^( Zl )dzi(l+ Zl ) 3 ^r(L, Zl ) X (11) 
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<7sis(£, zi,z s ) B(L, zi,z s ) x 



dL 

Id 



and the correspondence between L and 9 is given by Eq. (8) 



Note also that we have allowed, in these formulas, for a gen- 
eral dependence of the luminosity function, L, on zi, which 
we shall consider later in this paper. 



4.3 Dependence on parameters 

To illustrate the dependence of our observables (r and 
dr/dO) upon the parameters, we assume for a moment the 
following fiducial values: <j>* = 0.6 x 10" 2 ft 3 Mpc" 3 , a = -1, 
7 = 4 and cr» = 180km/s. For purposes of this illustration 
we have also assumed all sources to be at a fixed redshift, 
chosen to be z 3 = 1.3. 

The dependence of lensing statistics on the galaxy pa- 
rameters and various degeneracies between these parameters 
have been investigated extensively in the literature (see e.g. 
Kochanek 1993a, 1993b); here we present a brief overview. 
The variation of the total optical depth r around this fidu- 
cial model can easily be computed to be 

dlnr = 2.07dlnz s + l.OOdln0, +0.69dlna + 
4.16dlna* + 0.69 d In 7 - 0.61 din Q M + 
0.61 din w. (12) 

Perhaps not surprisingly, the strongest dependence is 
on the velocity dispersion, which strongly affects the lens- 
ing cross-section, as well as the luminosity function. <fi* en- 
ters linearly, and is degenerate with other factors, for exam- 
ple the magnification bias which is also a pure constant in 
the SIS case. Note, however, the much weaker dependence 
upon the cosmological parameters £Im and w. This rein- 
forces the notion, independent of observational uncertain- 
ties, that lensing constraints might most effectively be used 
to constrain galaxy profile and luminosity function parame- 
ters, in particular a*, rather than cosmological parameters. 
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Fig. 1 shows the dependence of r and (1/t){<1t /d6) on 
a,. As expected, r is a strongly increasing function of cr„, 
while (1 1 't){cLt / dff) favors higher angular splittings with in- 
creasing o t . As we shall describe, the fact that we only 
compare theory with observation for 9 > 0.3" (the angu- 
lar resolution of the survey) allows the likelihood function 
for angular splitting to be consistent with that for optical 
depth, which favors models with low cr*. 

We briefly comment on the dependence on other pa- 
rameters. Assuming 0* = const, only r depends on this 
quantity (we show in Sec. 9 that this is essentially true even 
if (j)* is redshift- dependent). Since r scales directly with 0*, 
the presence of other parameters implies that constraints on 
0, will be very weak. Furthermore, it is clear that, in the 
SIS case, r only depends on the combination a + 4/7 (this 
is slightly spoiled by the fact that the luminosity integral 
> 0) . We found that even trying to constrain 
this combination gives weak constraints - from either r or 
dr/d8 test. The only parameter that we are able to signifi- 
cantly constrain is a» . 



As mentioned in Sec. 2, we do not quote results from this test 
due to uncertainties regarding the redshift completeness. We 
do illustrate the constraints it gives in the GNFW case to 
demonstrate that including this result would not change our 
conclusions. 

The total likelihood we use is 

£tot = C T x Cdr/de (16) 

and it depends on cosmological parameters, as well as the 
Schechter function parameters </>», a, 7 and a». 

Equation (12) suggests that, for the SIS profile, by far 
the strongest dependence amongst the various lensing statis- 
tics is on the velocity dispersion <r* . We determine the like- 
lihood of a, by marginalizing over the other parameters: 

£(cr») = J £(<7», <f>*, at, 7) dcj)^ da dj. (17) 

where £ refers to any combination of the likelihood functions 
discussed above. 



5 THE LIKELIHOOD FUNCTION 

As we have mentioned there are several ways to use statistics 
of strong gravitational lensing. The total number of lenses - 
predicted vs. observed - is an obvious and most commonly 
used statistics which provides information about the inte- 
grated optical depth for lensing. The angular splitting and 
redshift distribution of lenses are the other statistical ob- 
servables, and in this work we use the former. We choose 
not to use the redshift distribution due to selection effects 
that are presumed to be significant in this test. Nevertheless, 
we have checked that the redshift distribution, if included, 
adds results consistent with the other two constraints. 

The probability of the total optical depth can be com- 
puted using the Poisson distribution (e.g. Kochanek 1993b) 



N x exp(-iV) 



(13) 



where x is the number of adopted lenses in the CLASS 
survey, N = 8958r is the number of galaxies predicted by 
the model, and r(0», a*, a, 7) is the computed optical depth 
given the Schechter function and cosmological parameters. 
This formula gives the correct likelihood for any value of r. 
Recall that our determination of r was based on the Marlow 
ct al. (2000) subsample redshift distribution. 

The likelihood for the angular distribution of galaxies 



-dT/d9 



1 dr 

n d6 



(14) 



where the product runs over the M lenses which we want 
to use for this test (recall, we use alternatively M — 4 or 
M = 9, and get virtually identical results for the two cases). 

Finally, the joint likelihood for the redshift and angular 
distribution of galaxies, which takes into account correla- 
tions between these two observables, is given by 



£d 2 r/dzd0 — Y]_ ~ 



dzd6 



(15) 



6 SIS PROFILE - RESULTS 

As mentioned above, the strong dependence of the opti- 
cal depth on the velocity dispersion cr* implies that we 
might hope to get an interesting constraint on a, despite 
the relatively poor lensing statistics and degeneracies be- 
tween lensing parameters. We marginalize over the other 
three relevant parameters, which we give top-hat (uniform) 
priors of <j>» G [0.5, 1.5] x 0.6 x 10~ 2 ft 3 Mpc~ 3 , 7 G [3.0,4.0], 
a G [—1.3,0.7]. These ranges are conservative, allowing the 
full spread of values reported in various recent measure- 
ments. We also made sure to use intervals that are sym- 
metric around the traditionally favored values, although it 
turns out that the exact choice of intervals affects the results 
very weakly. For example, the SDSS, from its commissioning 
data (Blanton et al. 2001), indicates that a = -1.20 ± 0.03 
and 0„,tot = (1.46±0.12) x 10~ 2 h 3 Mpc -3 , while the Two 
Degree Field survey from their preliminary sample of 45000 
galaxies (Cross et al. 2001) gives a = —1.09 ± 0.03 and 
0»,tot = (2.02 ±0.02) x 10~ 2 h 3 Mpc~ 3 ; both of these quote 
the total luminosity function. Kochanek et al. (2001), on 
the other hand, isolated early type galaxies from the K- 
band luminosity function, obtaining a — —0.92 ± 0.10 and 
<j>» = (0.45 ± 0.06) x 10~ 2 h 3 Mpc~ 3 . Finally, there are di- 
rect, independent constraints on the Faber- Jackson slope 7 
from the lens data; for example, Rusin et al. (2002) find 

7 = 3.44 ± 0.58 for early-type galaxies. 

Fig. 2 shows the 68% and 95% CL constraints from the 
r and drjdQ tests (top panels), as well as the constraints 
from the two tests combined (bottom panel). These con- 
straints correspond to solid curves in the three panels; for 
comparison, the dashed line in the first panel we show the 
effect of fixing </>„, a and 7 to their "fiducial" values, while 
the dashed lines in the top right and bottom panel indi- 
cate the effect of including the galaxies that are not iden- 
tified as ellipticals in the angular separation test. First of 
all, note that the two independent tests are in remarkable 
agreement, and that both constrain a* quite strongly. The 
r test gives 156km/s < a, < 226km/s (at the 95% CL), 
while the dr/d6 test gives 128km/s < a, < 272km/s (95% 
CL). Moreover, the dr/dO results are roughly independent of 
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a t (km/s) 

Figure 2. Constraints on the velocity dispersion a„ assuming the SIS lens profile and marginalizing over the other luminosity-function 
parameters. Top left: Constraint from the r-test (for comparison, the red-dashed line denotes the case when the galaxy parameters have 
been fixed to their fiducial values of tj>* = 0.6 X 10 — 2 h 3 Mpc — 3 , 7 = 4.0, a = —1.0). Top right: constraint from the dr / dd-test, assuming 4 
events that are due to ellipticals (solid line) and additional 5 events that are due to unidentified galaxies (dashed line). Bottom: the two 
tests combined; the two curves refer to the two subsamples used in the dr/dd test. Our baseline results, which we quote and to which 
the shaded confidence regions correspond, refer to solid curves in the three panels. 



the subsample of ellipticals we use, although the results are 
less tight in the baseline case when only the four "secure" 
ellipticals are used; see Fig. 2. The two tests combined give 
158 km/s < cr» < 220 km/s (95% CL). Therefore, the over- 
all favored value of cr* is actually smaller than the fiducial 
value of 225 km/s that has often been used to set constraints 
on cosmological parameters (Kochanek 1995, 1996, Falco, 
Kochanek & Munoz 1998, Waga & Miceli 1999, Cooray, 
Quashnock & Miller 1999) , and, not surprisingly, is in agree- 
ment with the value used in studies that tended to favor non- 
zero A (Cheng and Krauss, 1999, Chiba and Yoshii 1999). 
Note, however, that a* ~ 225 km/s has also been obtained 
using the direct observations of early-type lens galaxies (e.g. 
Koopmans & Treu 2002 get a, « (225 ± 15) km/s). Our re- 
sults disfavor this result as representing a fiducial value. 

We found that the constraint on <t„ is very weakly de- 
pendent on the exact value of intervals allowed for other pa- 
rameters. Furthermore, we find that independent constraints 
on other parameters of interest ((/>,, a and 7) are very weak, 
as expected from Eq. (12) and the fact that these param- 
eters are highly correlated (e.g. a and 7). Finally, we have 



checked that the dependence of these results on cosmology 
is extremely weak: for example, marginalizing over the plau- 
sible values of the matter density Qm € [0.15,0.40] (while 
maintaining the flatness condition) produces likelihoods that 
are only slightly broader. 



7 MODELING THE LENS: GNFW PROFILE 

There is a good evidence that galaxies have a cuspy inner 
profiles. The strongest argument comes from N-body sim- 
ulations, which argue for a profile p(r) oc r _/3 with (3 ~ 1 
(Navarro, Frenk & White 1996, 1997) or perhaps /3 ~ 1.5 
(Moore et al. 1999, Ghigna et al. 2001) - in either case, a rel- 
atively steep profile. Another argument in favor of strongly 
cusped central profiles is given by the absence of central im- 
ages in CLASS; assuming p(r) oc one obtains (3 > 1.8 
at 95% CL (Rusin & Ma 2001). Finally, direct modelling 
of the observed lenses favors steep inner cusps with profiles 
close to isothermal; p(r) oc r~ 2 (Munoz, Kochanek & Kee- 
ton 2001, Cohn et al. 2001, Treu & Koopmans 2002, Winn, 
Rusin & Kochanek 2002). These and other lines of evidence 
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suggest that the central profiles of lens galaxies are steep 
and that cores, if they exist, are tiny, with radius of a few 
tens or hundreds of parsecs at most. Such small cores would 
not affect the lensing observables appreciably (Hinshaw & 
Krauss 1987). 

To attempt to constrain the detailed profiles of elliptical 
galaxies we must move beyond the simple SIS model. In 
order to explore the dependence of lensing statistics on the 
details of the density profile, we adopt the generalized NFW 
profile described below. 



7.1 The GNFW profile 

The generalized NFW (GNFW) profile (Zhao 1996) is given 
by 



p(r) = 



I- 



3-0 



(18) 



where r s is the characteristic scale where the density profile 
shape can change. Because the integral of this density profile 
diverges at infinity, the mass of the halo is defined to be the 
mass contained within the radius r2oo at which the density 
is 200 times greater than the critical density of the universe 
at that redshift: 



' Air o 

W - Uxm = 200 ( — r 3 200 {z)p c (z) 



(19) 



The expression for the mass can further be written as 

f r 200( z ) 



= 4tt / 
Jo 



pr 2 dr = Airp s {z)rl (z) f{c{z)) 



M 

where 

/(c) = i ^(T+z)3-/3- 

and the concentration parameter is defined as 

r2oo(z) 



x dx 



c{z) 



r s {z) 

From Eqs. (19)-(22) it follows that 



(20) 



(21) 



(22) 



r 3 (z) 

Ps{z) 



3M 2I 



c(z) V 800717^(2) 

3 Pc{Z) f(c(z)Y 



1/3 



(23) 
(24) 



Thus, the generalized NFW profile is determined by the 
choice of the inner density slope j3 and the concentration 
c(z). Starting with these two parameters, one can compute 
p s (z) from Eq. (24) and then, given the mass of the halo, 
r s (z) from Eq. (23). Note that the GNFW profile for /3 = 2 
and the SIS profile are different for three reasons: 1) the 
GNFW profile parameters are explicitly redshift-dependent, 
2) the two profiles have different normalizations, and 3) the 
GNFW profile has a turnover at r = r s , while the SIS does 
not. 




Figure 3. The mean value of the concentration parameter as a 
function of the inner slope of the density profile f3 (solid line). The 
value at (3 = 1 and its l-cr uncertainty were obtained from N-body 
simulations. Concentration for other values of /3 was obtained by 
a simple recipe mentioned in the text, and adopting the same 
uncertainty in log 10 cq. 



7.2 The halo concentration 

The halo concentration factor c(z) is fortunately fairly well 
constrained due to recent results obtained using N-body sim- 
ulations (e.g. Bullock et al. 2001a, Wechsler et al. 2002). For 
a pure NFW profile, the concentration of the halos is well 
described by 



c{z) = 



CO 



+ 



M 



(25) 



with Co = 9 and M* = 1.5 x 10 13 Mq (the above papers actu- 
ally quote results for c vir = r v i T /r s with r v i r a virial radius, 
but the formula we quote accounts for the difference in defi- 
nition quite accurately) . The dependence on M is small and 
does not change the results much, while the dependence on 
redshift is important and fairly well-understood (Wechsler 
et al. 2002). It is also important to account for the variance 
in c which occurs not only because of uncertainties in halo 
modelling, but also because of the variance in halo proper- 
ties. We adopt an uncertainty in log 10 Co to be 0.14 (Bullock 
et al. 2001a, Wechsler et al. 2002). Therefore, when comput- 
ing the likelihood function we weight excursions around the 
middle value of Co by a gaussian factor with this standard 
deviation. 

Finally, we use the recipe from Li and Ostrikcr (2002) to 
compute co for a GNFW profile given c for a pure NFW: 
we assume that the ratio ri/ 2 /r2oo is independent of the 
density profile slope, where rj/ 2 is the defined as M(r < 
r \ji) — l/2M(r < 7"2oo). We retain the redshift and mass 
dependence of a GNFW profile as indicated in Eq. (25), as 
well as the same uncertainty in log 10 c . Figure 3 shows the 
mean value of the parameter Co and its standard deviation, 
both as a function of /3. 
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7.3 Cross-section for the GNFW profile 

Lensing by GNFW halos has been thoroughly explored by 
Wyithe et al. (2001) and Li & Ostriker (2002), and here 
we recapitulate the main results. The lens equation for a 
spherical symmetric lens is (Schneider, Ehlers & Falco 1993) 



(26) 



where j3 is the angular location of the source, 9 the angular 
location of the lens, and a the deflection angled D is and D s 
are the angular diameter distances between lens and source 
and observer and source respectively. Define £ and ff to be 
the position vectors in the lens and source planes respec- 
tively, and x = £/r 3 and y = (r//r s )(Di/D s ), where D t is 
the angular diameter distance to the lensing object. Then 
the surface mass density is given by 

-3+0 



roc 

s r s / (x 2 +z 2 )' 
Jo 



E(z) = 2p. 

and the mass by 

M(x) = 2-Kr 2 3 ( x'E(x')dx' 
Jo 



0/2 



((*W) 1/2 + l)) dz(27) 



(28) 



The deflection angle for a spherically symmetric source 



a(x) — 



4GM(i) 



The lens equation then becomes 



y = x- 
where 



M(x) 
47rp s rl 
4p s r s 



r 2 

Might 



4tvG DiD h 



(29) 
(30) 

(31) 
(32) 

(33) 



and ciight is the speed of light (to be distinguished from the 
concentration). Multiple images occur for x between ±x c , 
where x c is the solution of dy/dx — 0. Thus the cross section 
for the GNFW lens is 



(TGNFW 



■■ 7T [y(x c 



(34) 



7.4 GNFW optical depth for lensing 



The optical depth for the GNFW lens is completely specified 
by properties of the lens, /? and c(z), the locations of the lens 
and source, z t and z s , and the cosmological abundance of the 
lenses. As in the SIS case, we use the Schechter luminosity 
function to model the number density of galaxies, together 

§ The a and /3 used in this subsection are not to be confused with 
the Schechter function parameter a and the GNFW profile slope 
(3 used in the rest of the paper. Furthermore, note that ctgnfw 
is the cross-section, while a and cr* refer to the galaxy velocity 
dispersion. 



with the Faber-Jackson relation. The optical depth has the 
same form as in the SIS case: 



r(z s ) 



- 1: 

L 



dzi 



dL -rr{L, zi) a G NFw(zi,L) B(z t ,z 3 , L). (35) 
o aL 

In order to relate the optical depth to the parameters 
of a Schechter luminosity function, it is typical to define a 
one dimensional dispersion velocity of a GNFW profile in 
analogy to that defined for an SIS galaxy: 

2 GM 



2r 5 



(36) 



Combined with Eq. (19), this gives the mass as a func- 
tion of the dispersion velocity 



M G V WOnGpc 
This mass then determines r s 
1 V2 a 



r s (z) 



c(z) 10 H{z) 



(37) 



(38) 



which, together with p s {z), Eq. (24), specifies fx s , Eq. (32), 
which is necessary for the lensing equation. 

Finally, we need the magnification bias for the GNFW 
halos. For the source objects with the power-law flux distri- 
bution, as is the case with CLASS, this is given by (Li & 
Ostriker 2002) 

B = J^pAt 1 (39) 
where 



2x 



y(x c )y'(x )' 



(40) 



where prime denotes the derivative with respect to x and 
and xo is defined by 



y(x ) = 0. 



(41) 



Equation (40) has been adopted from Oguri et al. 
(2002); magnification bias defined this way agrees very well 
with ray-tracing simulations (C.-P, Ma, private communica- 
tion). The magnification bias for GNFW halos is very large, 
of order a few tens or hundreds. 



8 GNFW PROFILE - RESULTS 
8.1 Dependence on /3 

To compute lensing statistics using the GNFW profile, we 
need to supply c(z) and j3. Our main goal here is to deter- 
mine the inner density profile /3, which is a parameter of 
considerable interest and to which the lensing statistics are 
very sensitive. Therefore, we marginalize over the concen- 
tration normalization Co and parameters of the Schechter 
luminosity function. 

Figure 4 shows the total optical depth of a GNFW lens 
as a function of j3 for 1 < /3 < 2 and fiducial values of all 
other parameters. Also shown is the value predicted by the 
SIS model (horizontal dashed line), also with fiducial values 
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Figure 4. Optical depth vs. the value of the inner density slope, 
j3 (rising solid line) for the fiducial values of Schechter function 
parameters from Sec. 4.3. Also shown is the value predicted by 
SIS profile (horizontal dashed line). The optical depth inferred 
from CLASS is shown with the shaded region. 

of other parameters, as well as the optical depth actually 
measured by CLASS (shaded region). (As remarked before, 
there is no reason that the GNFW profile at (5 — 2 should 
match the SIS profile case.) From this figure it is clear that 
the optical depth is a strong function of (5. Note too that the 
lensing cross-section for > 2 is formally infinite, although 
it becomes finite if one considers configurations in which 
both images are detectable. While values (5 > 2 are allowed 
by our analysis, they are disfavored, and for computational 
reasons we only consider values /3 < 2. 

8.2 Parameter choices 

Our goal is to constrain the density profile /3. We therefore 
have to consider how to include a host of other parame- 
ters. We adopt the concentration function c(z) from N-body 
simulations, using Eq. (25) and choosing Co with a gaussian 
prior as discussed previously. We also need to marginalize 
over four luminosity function parameters ((/>,, cr*, a and 7). 
We choose the same ranges for a, 7 and (j>» as in the SIS case 
(see Sec. 6), plus a uniform prior a* G [150, 220] km/s, which 
is indicated by our SIS results. Remarkably, we find that in- 
teresting constraints on /3 are possible despite marginalizing 
over this large parameter space. As before, we assume the 
concordance cosmology (Qm = 1 — Ode = 0.3; w — —1). 

8.3 Constraints on f3 

The resulting constraints on the inner slope of the density 
profile are shown in Fig. 5. First, note that the total optical 
depth and angular separation tests (top panels) are in good 
agreement. The two tests together, when the likelihood func- 
tion is marginalized over other parameters, yield the con- 
straint 1.64 < p < 1.92 at the 68% CL and 1.50 < /3 < 2.00 
at the 95% CL (bottom left panel). As in the SIS case, these 



Figure 6. Same as Fig. 3, but overlapped with the 68% CL con- 
straint on /3 and co from the r-test of lensing statistics (filled 
circles). We use the total optical depth for the latter constraint, 
and retain the redshift and mass dependence of c(z) as in Eq. (25). 
Note that the region of overlap coincides roughly with the allowed 
range of these two parameters based on the likelihood function. 



results are insensitive to the exact ranges allowed for the 
luminosity function parameters. Moreover, we have checked 
that the angular separation test is insensitive to the choice 
of lens data, i.e. whether we use the 4 single deflectors con- 
firmed to be ellipticals, or all 9 single deflectors that are 
not identified as spirals. To be conservative, all results we 
quote correspond to the former choice and are represented 
by solid lines in Fig. 5 (for more on this choice, see Sec. 2). 
We also show the likelihood for the angular and redshift test 
combined using the three elliptical lenses with complete red- 
shift and angular separation information ((l/r)(d 2 r / 'dzid6); 
bottom right panel), which we did not use in the analysis 
due to uncertain systematic effects in the selection of lens 
redshifts. It is clear that the combined angular and redshift 
test is consistent with the other tests, and combining it with 
the r-test would further strengthen the final constraint on 
(5, as shown with the dotted curve in the bottom left panel. 

Although the favored slope is significantly steeper that 
the canonical NFW p oc r _1 profile, it is expected that the 
shallow NFW profiles seen in simulations become steeper 
due to baryonic infall (e.g. Kochanek & White 2001). The 
results of our analysis are in excellent agreement with such a 
scenario. Furthermore, these constraints are in good agree- 
ment with direct modelling of the observed lenses (Munoz, 
Kochanek & Keeton 2001, Cohn et al. 2001, Treu & Koop- 
mans 2002, Winn, Rusin & Kochanek 2002) which typically 
favors a steep, near-isothermal cusp. Finally, the results are 
insensitive to the exact values of cosmological parameters: 
for example, marginalizing over the plausible values of the 
matter density Qm £ [0.15,0.40] produces negligible in- 
crease of the width of our contours. 

We also can constrain the GNFW concentration param- 
eter c(z) and the density profile slope f3 jointly. In Fig. 6 we 
display the N-body determination of the concentration pa- 
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Figure 5. Constraints on the inner density slope j3, marginalized over all other relevant parameters. Shown are the likelihood functions 
for r and (l/r)(d,T /d9) (top panels), as well as for the two combined (lower left panel). The lower right panel shows the likelihood for 
the angular and redshift test combined ((l/r)(d 2 r/dzidB)), which we did not use in the analysis due to the uncertain redshift selection 
function, but show for illustration that it is consistent with the other tests. For the likelihoods using the angular separation test, we 
show the results assuming 4 single deflectors confirmed to be ellipticals (solid lines), and, alternatively, all 9 single deflectors that are not 
identified as spirals (dashed lines). Note that the solid and dashed line in the combined likelihood test esentially overlap. Dotted line in 
the combined likelihood test shows the likelihood when r and (l/r)(d 2 r / dzidO) tests are combined. 



rameter as a function of /3, and overlay this with our lensing 
constraint on cq vs /3, using the r-test. For any given /3, we 
allow Co to be a free parameter, and retain the redshift and 
mass dependence of c(z) as in Eq. (25). Not surprisingly, 
the allowed value of /3 reported above coincides with the 
overlap region between the N-body result and our lensing 
constraint. Note, however, that lensing imposes constraints 
on the concentration that are independent of N-body results. 
In particular, if the galaxies indeed have pure NFW (/3 = 1) 
profile, lensing statistics implies that the concentration pa- 
rameter Co has to be greater than 15, which is in conflict 
with the results of N-body simulations. 



9 REDSHIFT-DEPENDENCE OF THE 
LUMINOSITY FUNCTION? 

We mentioned previously that one of the great difficulties 
with using gravitational lensing statistics as a probe is that 
the parameters that describe the abundance of galaxies can 
depend on redshift. (In the GNFW case, the concentration 



parameter c(z) is allowed to vary with redshift, as predicted 
by numerical simulations. ) To make progress, essentially 
all authors in the past who wanted to use lensing statistics 
assumed that these functions were redshift-independent. In 
particular, one expects that the number density <^>* and the 
characteristic velocity dispersion a* may be strongly depen- 
dent on redshift due to galaxy accretion and mergers. ^. 

Direct constraints on the redshift dependence of the lu- 
minosity and abundance of galaxies are still crude, made 
difficult by poor statistics and a variety of systematic ef- 
fects. Even rough agreement between various surveys has 
not been achieved. For example, while the Canada-France 
Redshift Survey (CFRS; Lilly et al. 1995), the CNOC2 sur- 

^ This situation is reminiscent of that in the analysis of galaxy 
surveys, where one needs to know the galaxy-to-mass bias in order 
to obtain the distribution of matter from the observed distribu- 
tion of galaxies. In the past most authors assumed the bias to be 
constant, while it is widely suspected that it depends on scale, 
redshift, and galaxy type. 
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vey (Lin et al. 1999), and the CADIS survey (Fried et al. 
2001) all observe an increase of 0* for early-type galaxies 
between redshifts of zero and z ~ 1, the Autofib survey (El- 
lis et al. 1996) and the numerical simulations by Nagamine et 
al.(2001) conclude just the opposite. It is clear that getting 
the redshift dependence of number densities and characteris- 
tic velocities per spectral type and their various covariances 
will take some time. Keeton (2002) has argued that a vari- 
ation in with 2 can cancel out much of the cosmological 
sensitivity of lensing statistics. However, we note that this 
variation alone is probably unrealistic. At the same time, 
mergers and accretion will be expected to cause a variation 
a, , which will have the opposite effect of a variation in <j>„ 
on lensing statistics, and indeed may overwhelm it. To ac- 
curately account for evolution, it is probably best to match 
onto N-body simulations of the galaxy mass function, which 
in fact suggest that the number density of galaxies with a 
specific value of cr* is relatively constant with z (i.e. Bullock 
et al. 2001b). 

To estimate the maximal possible effect of evolution (as- 
suming an evolution in the galaxy number density only), we 
used an SIS profile, which simplifies calculations. If one then 
considers a number density dependence of galaxies as sug- 
gested by Lin et al. (1999) 

Mz) = <M0)10°' 4Pz (42) 

one can estimate how the results would change for non-zero 
values of P. For the r-test, the change is as expected: for 
example, for P = 1 the number density increases by ~60% 
(assuming the average lens redshift is ~ 0.5), which corre- 
sponds to the decrease in the favored cr* by ~10% in order 
to preserve agreement with the measured r, cf. Eq. (12). For 
the {1/t)((Lt /dff) test, the redshift-dependence of 4>* largely 
cancels out in the numerator and denominator of this quan- 
tity. Therefore, as expected, the total optical depth is more 
sensitive to the redshift dependence of 4>, , while the angular 
distribution of lenses is not. Again, we expect that the ac- 
tual impact of evolution will be much less severe than that 
discussed above, because mergers and accretion will tend to 
produce a variation in u* with z that will cancel the effect 
of the variation in <^>*. 



10 CONCLUSIONS 

The use of strong gravitational lensing statistics in order to 
probe cosmology has a long history. Nevertheless, the dom- 
inant uncertainty in the predictions of lensing statistics has 
to do with estimates of galaxy parameters, not cosmological 
ones. Because of the recent revolutions in observational cos- 
mology that have allowed us to pin down the basic cosmolog- 
ical parameters with relatively good accuracy, gravitational 
lensing statistics now provide us a new opportunity to probe 
the structure of galaxies and the trends of galaxy evolution. 
Our results represent a first step in this regard. Neverthe- 
less, it is quite remarkable that, in spite of the paucity of 
lensing statistics at this time, we obtain non-trivial limits 
on galaxy properties. It is also significant that these limits 
are largely independent of cosmological uncertainties. Since 
we are primarily interested here in constraining observa- 
tional galaxy parameters we used the Schechter luminos- 
ity function, which gives the number density of galaxies in 



terms of luminosity, rather than the mass function, which is 
more relevant for more massive halos associated with clus- 
ters (M > 1O 13 M ). 

Assuming the SIS density profile, we find that the 
mean velocity dispersion for elliptical galaxies is small, with 
168km/s < a, < 200km/s at 68% CL, consistent with a 
number of earlier estimates used in lensing analyses (e.g. 
Chiba & Yoshii 1999, Cheng & Krauss 2001), but signifi- 
cantly smaller than the "canonical" value of 225 km/s often 
quoted in the literature. Perhaps more significantly, assum- 
ing the generalized NFW density profile with inner slope 
1 < /3 < 2, we constrain (5 to be in the range 1.64 < /3 < 1.92 
at 68% CL. This is definitely inconsistent with the j3 — 1 
slope advocated by N-body simulations for the dark matter 
halos profiles; Fig. 6 shows that a profile with (5 = 1 could 
produce the observed lensing statistics only with an unrea- 
sonably high concentration (co > 15). At the same time, it 
is a well-known fact that N-body simulations do not include 
additional physics, e.g. the baryonic infall, that makes the 
inner profiles of halos and galaxies steeper. Consequently, 
our result for the density slope is in good agreement with 
the expectations, as well as with similar analyses (e.g. Kee- 
ton 2001, Kochanek & White 2001, Rusin & Ma 2001) or 
direct modelling (Cohn et al. 2001, Treu & Koopmans 2002, 
Winn, Rusin & Kochanek 2002). The lack of high-separation 
events (>3") in JVAS/CLASS has been used to suggest that 
there are two populations of halos in the universe (Keeton 
1998, Li & Ostriker 2002): small-mass galaxy-size halos with 
possibly steep density profiles (/3 ~ 2), and large- mass halos 
with shallow density profiles (/3 ~ 1). The former correspond 
to the elliptical galaxies we are interested in here, and our 
results confirm that a steep slope seems to be required to 
explain these events. 

There remain some issues that require further explo- 
ration. In particular, the total optical depth produces a like- 
lihood function that tends to suggest a slope that is less steep 
than that favored by exploring the redshift dependence and 
angular splitting of lensing events. This may be an artifact 
of our limited statistics, but it could also signal the need 
to consider a more complicated, perhaps two component, 
galaxy distribution in order to consistently model lensing 
events. 

Needless to say, the most significant factor not explic- 
itly taken into account here is a possible redshift evolution 
of galaxy number density and velocity dispersion. Very lit- 
tle is known observationally about the evolution of these 
quantities beyond z ~ 0.3. As we have discussed in Sec. 9, 
there are reasons to believe that the effects of evolution will 
not significantly alter the allowed parameter ranges we have 
determined here. 

Finally, we note that there is great potential to im- 
prove these constraints as better statistics are obtained us- 
ing current and future observational efforts. In particular, 
the DEEP2 redshift survey (Davis et al. 2002) will provide 
a velocity function for galaxies at redshift z ~ 1, which will 
allow one to explore the evolution of galaxy parameters with 
redshift with a much higher sensitivity than currently avail- 
able. 
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